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Abstract 

We derive the radiation power from a single Josephson junction (JJ) and from layered superconductors 
in the flux-flow regime. For the JJ case we formulate the boundary conditions for the electric and magnetic 
fields at the edges of the superconducting leads using the Maxwell equations in the dielectric media and 
find dynamic boundary conditions for the phase difference in JJ which account for the radiation. We derive 
the fraction of the power fed into JJ transformed into the radiation. In a finite-length JJ this fraction is deter- 
mined by the dissipation inside JJ and it tends to unity as dissipation vanishes independently of mismatch 
of the junction and dielectric media impedances. We formulate also the dynamic boundary conditions for 
the phase difference in intrinsic Josephson junctions in highly anisotropic layered superconductors of the 
Bi2Sr2CaCu208 type at the boundary with free space. Using these boundary conditions, we solve equa- 
tions for the phase difference in the linear regime of Josephson oscillations for rectangular and triangular 
lattices of Josephson vortices. In the case of rectangular lattice for crystals with the thickness along the 
c-axis much larger than the radiation wavelength we estimate the radiation power per unit length in the 
direction of magnetic field at the frequency 1 THz as fiW/cm for Tl2Ba2CaCu208 and 0.04 N 
/iW/cm for Bi2Sr2CaCu208. For crystals with thickness smaller than the radiation wavelength we found 
that the radiation power in the resonance is independent on number of layers and can be estimated at 1 THz 
as 0.5 W/cm (Tl2Ba2CaCu208) and 24 mW/cm (Bi2Sr2CaCu208). For rectangular lattice, due to super- 
radiation regime, up to half of power fed into the crystal may be converted into the radiation. In the case 
of triangular or random lattice in the direction perpendicular to the layers the fraction of power converted 
into the radiation depends on the dissipation rate and is much lower than for rectangular lattice in the case 
of high- temperature superconductors with nodes in the gap. 
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INTRODUCTION 

In 1962 Josephson 111 predicted electromagnetic radiation from superconducting tunneling 
junction arguing that in the presence of voltage V across the junction the phase difference 
Lf changes with time as dip/dt = 2eV/h, while the tunneling current density J changes as 
J = JcSm(2eVt/h). Here Jc is the critical superconducting current via the junction. Thus the 
photons with the Josephson frequency u = 2eV/h may be emitted from Josephson junction (JJ). 
Such a radiation from J J into the waveguide with the power Vrad ~ 10"^^ W has been detected 
by Dmitrenko et a/.|Q] and by Langenberg et al. isl soon after the Josephson prediction. The 
junction in both experiments was subject to bias current and dc magnetic field applied parallel to 
the junction. The dc magnetic field induces Josephson vortices and in the presence of the trans- 
port current across the junction they move to the JJ edge causing oscillations of the magnetic 
and electric fields inside JJ and in superconducting leads around JJ. More precisely, flux flow in- 
duces electromagnetic (Swihart) waves inside JJ, and they are partially transmitted outside when 
hit JJ edges lO]. From measurements of the current- voltage (I-V) characteristics, Langenberg et 
al. Isl found the power V = IV fed into the junction and the fraction converted into radiation, 
Q = 'Pva.d/'P ~ 10"^. Theoretically Q was estimated as a ratio of the impedances of JJ modeled 
as a strip line and the waveguide ||3|,|5|]: 

where Zq is the waveguide impedance, Zs is the strip line impedance in the limit k^^w <^ 1 
and Zg <^ Zq. Here A is the London penetration length of the superconducting leads, d is the 
thickness of the insulating layer, 6i is its dielectric constant, = cu/c and w is the junction width. 
Eq. (UJ) gives ^ 10"^ for the studied junction in agreement with experimental value. Thus low 
radiation power from JJ was attributed to the mismatch between the impedances of the junction 
and the waveguide. To the best of our knowledge, no deeper theoretical treatment of radiation was 
made after that and the concept of the impedance mismatch for a single JJ and later for layered 
superconductors with intrinsic JJ was accepted by the community. 

As radiation from the vortex flow in a single junction was found to be quite low, a natural 
idea was proposed to use multiple lock-injunctions. Then, in the super-radiation regime, radiation 
power may be enhanced by the factor A^^, where N is the number of synchronized junctions inside 
the space of the radiation wavelength. Intensive theoretical and experimental study was done in 
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this direction, see e.g. Ref. However, an effective way to synchronize many 

junctions was not found so far. There are two problems to overcome for synchronization. The first 
one is the technologically inevitable variation of parameters from junction to junction, mainly 
parameter, which affects operating frequency of the junction at given current. The second problem 
is that the coherent locked-in flow of the Josephson vortices in different junctions is unstable in a 
wide range of parameters. 

The discovery of layered high-temperature superconductivity added new breath into this activ- 
ity. It was recognized that the Bi- and Tl-based cuprate superconductors with weakly coupled su- 
perconducting Cu02 layers exhibit the same static and dynamical Josephson properties as artificial 
tunneling junctions. In another words, the crystals Bi2Sr2CaCu208 (BSCCO) and Tl2Ba2CaCu208 
(TBCCO) represent a stack of many intrinsic Josephson junctions on the atomic scale lisi . The 
first indication of the intrinsic Josephson effect was observation of switching of individual junc- 
tions to the resistive state in the I-V characteristics 13]. After that many Josephson effects have 
been found in these systems including the Fraunhofer patterns in the dependence of the critical 
current on the dc magnetic field applied parallel to the layers lisl . Josephson plasma resonance 
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ular to the layers lll7Lll8ll . Shapiro steps in the I-V characteristics induced by external microwave 



Il6i strongly affe cted by pancake vortices in the presence of the magnetic field applied perpendic- 
e lave 

radiation ligi lioi. and Fiske resonances Il21 . 

One can anticipate much smaller variations of intrinsic-junction parameters in comparison with 
the artificially-fabricated junction arrays. In addition, the intrinsic junctions are much closer to 
each other and one can anticipate much stronger coupling between them. What is more, there 
are many junctions on the scale of the radiation wavelength, and so they will super-radiate when 
synchronized. Due to these advantages, the moving Josephson vortex lattice in layered super- 
conductors was proposed as a source of monochromatic tunable continuous electromagnetic ra- 



diation in the terahertz frequency range 



2811 . Experimentally, radiation from the 



high-temperature layered superconductor BSCCO at relatively low frequencies, 7-16 GHz, was 
detected by Hechtfischer ei al. Some indirect evidence of radiation at higher frequencies 

has been recently reported by Kadowaki et al. jsf] ] and Batov et al. Isil reported radiation at 
the frequency 0.5 THz with power 1 pW from BSCCO mesa consisting 100 junctions in zero dc 
magnetic field. 

Therefore, it is interesting to estimate theoretically the possible radiation power generated by a 
moving vortex lattice and find optimal conditions for generation. For that we need to understand 
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the mechanism of conversion of the electromagnetic field associated with the flux flow of the 
Josephson vortices into the electromagnetic waves outside of the Josephson junction (JJ) and find 
limitations imposed by this mechanism in Josephson tunneling structures. A natural first step is 
to understand at the microscopic level the conversion mechanism in a single JJ. Then such an 
approach may be extended to layered structures. 

In the first part we consider radiation from a single JJ. Then the method to treat radiation is ex- 
tended to intrinsic junctions in layered superconductors. We will show that super-radiation regime 
is inherent to moving rectangular vortex lattice in such crystals. We will discuss the consequences 
of this regime for the radiation power and I-V characteristics. 



RADIATION FROM A SINGLE JOSEPHSON JUNCTION 

As we discussed in the Introduction, a common way to evaluate the radiation out of a single 
Josephson junction is to take product of the total power supplied to the junction, V = IV and 
the impedance mismatch coefficient Qz, Eq. (HJl. A weak point of this impedance approach for 
real, finite-length JJ is ignoring of multiple reflections inside the junction, i.e., assumption that 
the propagating electromagnetic wave has only one attempt to escape the JJ, decaying before 
reflected wave reaches another edge. In fact, at low dissipation rate (low temperature and low 
dielectric losses inside the insulating layer) reflections lead to the formation of almost standing 
Swihart waves inside JJ. In this case Q strongly depends on the dissipation and approaches unity 
as dissipation vanishes. Then the question turns out to be what are limitations on V and VraA 
rather than on Q. 

On the other hand, the standard analysis of transport properties of finite-size JJs uses sine- 



jquation and zei 



Gordon equation and zero-derivative boundary conditions /or the oscillating part of the phase at 
the edges 



At the JJ edges vanishing oscillating magnetic field B = d^pjdx leads 
to zero Poynting vector. Thus for such boundary conditions the outside radiation is absent and all 
supplied power dissipates inside JJ. 

In the following we reconsider this problem, discuss the power conversion mechanism and 
derive the radiation power from a single JJ into free space in the case when w is much larger than 
the wavelength of outcoming electromagnetic wave. This is not realistic limit, but it is a simplest 
one which demonstrates the method to treat radiation and which may be directly extended to 
the case of smaller w. Our rigorous approach is based on solution of the Maxwell equations 
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FIG. 1: A finite-size Josephson junction opened into free space at both edges. Ellipses illustrate the moving 
vortex lattice. Curved lines show the screening currents inside the superconducting leads, arrows show 
radiation from area 2A around the junction. The applied dc magnetic field Hq is along the z-axis. 

inside the superconducting leads and in outside space, which allowed us to formulate accurate 
dynamic boundary conditions for the oscillating phase inside JJ. In the linear regime of Josephson 
oscillations we obtained analytical results for "Prad and Q using the perturbation theory. In this 
regime we found that Q oc Qz^ ■, where M is the number of reflections before Swihart wave 
decays inside JJ. At low temperatures in JJ made out of gapped superconductors with perfect 
insulating layer M may be large compensating small Qz- Our approach also opens the way for 
numerical calculations in general case, when linear approach is invalid. 

The Josephson junctions with low level of the dissipation become available now due to the 
perspective to use them as a qubits for quantum computing 113 5h . The radiation should be stronger 
in such junctions in comparison with previously studied ones. This gives additional motivation to 
reconsider the theoretical background of the radiation from JJ. 

To find outside radiation due to the Josephson oscillations, one has to match oscillating fields 
inside the junction and in the superconducting leads with the wave solution outside the junction. 
Expressing the oscillating fields inside the junctions via the phase difference, we derive dynamic 
boundary conditions for the phase difference and relate the Poynting vector of radiation to the 
phase difference at junction edge. The final step is the solution of the equations for the phase 
difference accounting for the dynamic boundary conditions and derivation of the Poynting vector. 



Equation for the phase difference 

We consider a JJ with the length / ^ A located at — / < z < and bounded by dielec- 
tric with dielectric constant e^, see Fig. [B Strength of the coupling in the junction is char- 



5 



acterized by the Josephson current density Jc and related parameters: the Josephson length, 
Aj = A/c$o/(8vr2AJc), and plasma frequency, ujp = ^/Sir'^dcJc/iei^o), where A = 2X + d. 
We consider the simplest situation when the JJ width along z direction, w, is larger than both Aj 
and wavelength of outcoming electromagnetic wave. We will consider straight Josephson vortices 
along the external magnetic field. In this approach the problem becomes one-dimensional for the 
phase difference and two-dimensional for the electric and magnetic fields in the outer space as 
both of them become z-independent. We consider junction in resistive state and assume that the 
junction phase, (/?, oscillates with the Josephson frequency u generating oscillating electric (E^ 
and Ey) and magnetic (B^) fields inside the junction and in the superconducting leads. Our task is 
to find spatial distribution of these fields and match them with outside fields to find equation and 
boundary conditions for 

We first derive equation for the oscillating magnetic field inside the superconducting leads. We 
use complex representation for the oscillating fields and phase, e.g.. 



(f{x, t) = {(f{x, t))t + ^ Re[v9^^(x)e" 



iujt] 



Here {■ ■ .)t means time average. The phase gradient is connected by the following relation with 
the magnetic field inside the junction and supercurrents flowing along the junction at the opposite 
sides 

V.^ = [J^^ - J^-] - -^B.. (2) 

From Maxwell equations, material equation for current inside superconductor, J = 
(c/47rA^) [(<l>o/27r) V0 - A] + cTgE, London relation for the electric field, E = - (AnX'^iuj /c) J, 
and Eq. Q, we derive the following equation for the oscillating magnetic field inside the leads 

(-/ < X < 0) at d < A 



XI 27rA2 dx 

Aj^ = A"^ - Esk^ + A-Kik^a^/ c, k^ = u/ c. (4) 



(V, + V,J B-,{uj) - — — = ^^^5{y), (3) 



We ignore small contribution to the magnetic field from the dc current flowing via the junction. 
The ac electric field inside the superconducting leads is Ey{uj, x, y) = iX^^k^V z{ijJ , x, y). To 
obtain the total electric field in the junction area, one has to account also for the field inside the 
dielectric layer. This gives 

Ey(uj, X, y) = ~^-^^^if^{x)5{y) +iXlk^'S/ .^Bz{uj, x, y). (5) 



The boundary condition for ipu}{x) follows from the boundary condition for the electric field 

-iuj[edE.,{+Q, y) - eE,{-0, y)] = 47rJ,(-0, y). 

As 0, y) = —'ico'(47rA^/c^) Jj;(—0, y), we obtain 

47r(l - eklX^)J,{-0, y) = -zo;£rfE,(+0, y). (6) 

As the electric field i?a;(+0,?/) is continuous at y = 0, this means that Jx{0,y) also must be 
continuous and from Eq. Q in the limit o? <^ A we obtain 

V^if^O) = V,.^^(-/) = 0. (7) 

We can represent solution of Eq. ^ near the edge x = as 

Bz{x,y) = Bb{x,y) - ^^^^^^ J Wx'^u;ix')dx' x (8) 
[A^o(c_) - A^o(c+)], c± = ^/{x±x')^ + yyX^, 

where Kq (z) is the modified Bessel function and -Bfe(x, y) is the solution of the homogeneous 
equation 

{Vl + Vl)B,-\z'B, = 0, (9) 

with the boundary condition Bb{—0,y) = Bz{+0,y). As a function of x, Bb{x,y) decays at 
distance ~ A from the boundary. For its Fourier transform along the y direction, we obtain [k^ = 

Bb{uj,x, ky) = B^{+0, ky) exp(-K^.|x|). 

For static case, distribution of the magnetic field near the edge of JJ has been recently derived by 
A. Gurevich (unpublished). 

Using Maxwell equation V = {iuJSi/c) Ey — (47r/c) Jy, Josephson relation Ey{x, 0) = 
—iuj^Qtp^^iyx) I {27icd), and Eq. we find equation for the oscillating phase 

— + at— <^Ux) + -i / dxVx Ko — + 

^pj ttA L V / 

_ cVxBb{x,0) 

^x'^uj{x ) = S^{X), (10) 



V 

X \ x' 



where s^(x) is the complex amplitude of the oscillating Josephson current, sm[if{x,t)] = 
Re[suj{x) exp{—iiut)], and at = ujpat^o/ {^ixJccd) is damping due to the tunneling quasiparticle 
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conductivity. Similar nonlocal equation has been derived by Gurevich ll36ll . In the case A ^ Aj 
and |x| ^ A non-locality is not essential and Eq. (flOb can be reduced to usual local approximate 
equation 



-at- 



UJr 



(11) 



where ag = 27r X'^lu pa g/c^ is the dissipation due to quasiparticles inside superconductor and we 



used expansion A^/A 



27r«A^A;a;0-5/c neglecting /c^ term. 



Boundary conditions and the Poynting vector 

Near the boundary situation is more complicated because, in addition to smoothly changing 
part, the phase has component decaying at distances of order A from the boundary. This extra 
phase is smaller than the smooth phase by the factor ~ A/Aj but it has comparable derivative. 
Our purpose is to derive accurate boundary condition for the smooth phase, obeying Eq. 

dTTTl . For this we integrate Eq. dlOb from intermediate distance —Xi with A <^ Xj ^ A j up to 
the boundary a; = 0. Neglecting small terms proportional to Xi and Bb(—Xi, 0) [for -82(0, y) <^ 
$o/(47rA^)], using local approximation at x = —Xi, we obtain the boundary condition for the 
smooth phase in a very simple form 

V.<^^(0) ^ V,^^{-Xi) = -^^B,{iu,r = 0). (12) 

I'd 

This equation can be compared with condition Q for the total phase. Eq. (fT^ allows us reduce 
the problem to solution of local equation (fTTTl for smooth phase with modified boundary conditions 
and avoid solving exact integral equation for the total phase, Eqs. dTJ and (flOb . 

For that we need to express -^^(r = 0) via ip^. Relation between the electric and magnetic field 
at the boundary is determined by properties of outside media. 



Outside dielectric media 



We consider first outside dielectric media with dielectric constant sa- In the situation wk^^ ^ 1, 
in the straight-vortices approach, outside fields are ^-independent. In addition, if we assume that 
the dielectric media is infinite in y direction and the thickness of the leads is much larger than A, 
we can use Fourier transform in this direction. In this case the Fourier components of fields with 



8 



\ky\ < y/£dkuj propagate in the media while the field components with \ky\ > ^Jelk^^ decay. In 
particular, for Ey{uj, x, ky) we have 



,^ (Ey{u, 0, ky) exp[i^edkl - k'pign{uj)x], for \ky\ < ^/e^k^, 

Ey{uj,x,ky) = < ^ (13) 

t Ey{uj,0,ky)exp[-^yk^ - Edklx], for \ky\ > y/e^\k J . 

Other field components, E^ and B^, can be expressed via Ey(uj, 0, ky). First, Ex{uj, x, ky) can be 

obtained from Eq. (fOl) and Maxwell equation V E = and then Bz{uj, x, ky) can be obtained 

using the Maxwell equation (V x E)^ = ik^^Bz which gives 

£d\ku 



Bz{uj,x,ky) = —====Ey{uj,0,ky)exp[iJedkl-k'^sign{uj)x], for \ky\ < ^/e2k^, 
V ^dk^ — ky V 

Bz{uj,x,ky) = ^^'^^^^^ Ey{u,0,ky)exp[-Jk^-edklx], for \ky\ > y^k^, (14) 
yky— £dk^ V 

In particular, this gives relation between fields at the boundary, which we will use to formulate the 

boundary conditions for the phase 

5,(0, ky) = C{UJ, ky)Ey{0, ky) , (15) 

(^(^ ^ ) _ / Ikjsd/ ^/edkl^^Jor \ky\ < ^/e2\kj, 
^ X-ik^Ed/ ^Jkl - Edk'i, for \ky\> ^Je2\k^ \ . 

Note again that the term C{uj, ky) for \ky\ < ^JWdk^^ originates from outcoming electromagnetic 
wave (radiation), while the term k^ for \ky\ > ^Jedk^ is due to the wave decaying at dis- 



tance ~ (/c^ — e^/;;^)"^/^ from the lead boundaries. The latter term does not carry energy out 
of the junction. For completeness, we also present this important relation in the frequency-space 
representation 

B^uj, 0, y) = y dy'\J{yj, y - y')Ey{u;, 0, y'), (16) 
U{uj,y) = ^ [\kJJo{^/e2k^y) + ik^No{y/e2k^y)] . 



2 

where Joiz) and Nq{z) are the Bessel functions, and in the time-space representation 

B,{t, 0, y) = v/^l- / dt'dy' '^^ ';^^^^~^'^'~^y~y'^'^ Ey{t', 0, y') (17) 

where 6(x) = 1 if x > and if x < and Cd = c/ ^/Sd- 

Now we will relate boundary fields with the phase. From Eqs. © and © follows the relation 

Ey{Q, ky) = tXlk^K,B,{0, ky)- (18) 



2n 
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This relation and Eq. (fTSb allow us to express the boundary fields via the phase distribution 

$0 iCku. 



2n 1 - i(\lk^K, 



■X 



(19) 



Typically, C-^^^w'^a; ~ Xtu/c <^ 1. Also V2:(/?a;(a;) ~ ^uj{x)/Xj and the integral term in Eq. (fT9b is 
smaller than v^^jIO) by the parameter A/Aj. Therefore, we obtain 



$0 ^0 

Ey{0, ky) ^ -—ik^^^{Q), B,{0, ky) ^-—iCk^ip^iO). 



(20) 



27r \ 27r 

Within this approximation, the magnetic field at the junction edge which determines the boundary 
condition (fT2b is 



5,(r = 0) ^ -5^zA;^(^^(0)Z(aj), Z(c^) 



(21) 



^ dky £d\ku. 

-^Ci^, ky) ^ — - 



.Edk^ 71 /d 
I m 



vr \fe'd\kj\ 

Finally, we obtain the boundary conditions for smooth oscillating phase at both edges, a; = 0, — /, 
in a finite-length JJ: 



Va;V5^(x) ^ ±H\k^Z{iJ)if^[x),ioxx = 0,-1. 

Radiation outside JJ to the right is given by the Poynting vector, Pj-ad 
ic/A7i)Jdy{Ey{0,y,t)B,iO,y,t))t, 



(22) 



(23) 



We also present the phase equation and boundary conditions in the time representation. Introduc- 
ing the dimensionless variables r = tUpt and u = x/\j and using Eq. (fTTT) . we write the equation 
for the total smooth phase ^^(r, u), which includes both oscillating and non-oscillating components 



52 g 



ip + sirup = 0. 



(24) 



Using Eqs. ([T2l) . (1211) . and (l22b and adding the boundary condition for the static phase, we obtain 
the dynamic boundary condition: 



VuVir, 0) 



rfr'/C(r - t' 



Or' 



(25) 



/C(r) 



dzJo{z). 
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FIG. 2: The Josephson junctions at a; < 0, dielectric at < x < L ^ k~ and superconducting screen 
with the London penetration length As and the thickness D ^ \s. Such a screen effectively reflects the 
Swihart wave back into the junction. It introduces additional interaction of vortices in different junctions in 
the case of multiple pai^allel junctions. 

Here, ^ = nc/ {cUpdy/Sd), Jo{z) is the Bessel function, b = Ah\\jHq/^q with Hq || z being 
the applied dc magnetic field, h <^ Aj/A. For the edge x = — / we need to change sign of /C. 
Due to the non-analytical cj-dependence of Z{uj), the kernel /C(r) is irregular at ^ and then 
r ^ 0. Singular frequency dependence of Z{uj) is due to retardation caused by electromagnetic 
wave propagation inside dielectric, Eq. (fTTl) . Thus boundary conditions also exhibit retardation 
effect. 



Outside superconducting screen 

Let us consider now a junction with superconducting screen at the right side separated from the 
layered crystal by the dielectric with the thickness L small in comparison with the wave length 
of the electromagnetic wave in this dielectric, see Fig. |2| Such a screen prevents the radiation 
from the right side increasing it from other side. In the case of many parallel junctions the screen 
enhances their interaction as was suggested by Ivanchenko flli . 

We consider the case k^L <^ 1 and L <^ Xg. Fields inside the superconducting screen obey 
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Eq. The solution is 

1 f~^°° ik 
E^ico, X, y) = -— / dkyG{ky)^e-^^^^^''yy, (26) 

2 



1 /■+°° 



where A^^^^, given by Eq. ©, characterizes the superconducting screen and Xg is the London pen- 
etration length of the screen. From this solution we obtain 

B,{uj, r = 0) = -%k^ip^{Q)Z,{uj), (27) 

ZTT 



i , /2A 



^ Tik,„\l.,. V L 
This leads to the boundary condition 



- W..(0), ^ ~ —r^ In — . (28) 



du TiXl \ L 

Note that 77 ^ 1. This results in the condition that V5a;(0) is very close to zero for screened edge of 
the junction, while at the open edge the condition is that space derivative is close to zero. There is 
flow of energy from closed edge proportional to the small parameter AixX^^aq/c due to quasiparticle 
dissipation inside the screen. 

Solution for the phase difference, I-V characteristics, radiation and dissipation power 

Now we solve analytically Eq. (l24l) using the perturbation theory with respect to the Josephson 
current in the limit |6 — cj|6 ^ 1. Taking solution as 

ip{T,u) = UJT -hu + e{T,u), 6'(r,M)<l, (29) 

and expanding sm[^p{T, u)], we see that 9^{u) = 9{0j, u) obeys reduced versions of equation (fTTl) 
with = -e^^/i, 

[Wl + Cj^ + iCj{at - a,V,'J] = -e^'^'/t. (30) 
and boundary conditions (|22|l 

^ = ±iCu/3ip, for u = 0, -I (31) 
du 

12 



with 



13 = l3o\ \uj\ 111' 



TT 



Solution for 9{uj, u) given by, 



— + aie*^"'' + 026 



(32) 



_ ^2 _j_ {q,^[^ 

describes the moving vortex lattice (the first term) and reflected Swihart waves propagating to 



the right and left. Here = at + Ogb"^ and = a/ {(^"^ + iatuj) / (1 — iaqOj) ^ Cj + ia/2 with 
a = at + aqUj'^. For a JJ with both open edges finding ai and 02 from the boundary conditions (|3T1) 
we obtain the phase distribution 

exp (ibu) /i 



9Ju) 



^2 _ ^2 _|_ ^q,^lD 



(b-Pu) 


cos[p^{l+u)]-if3 sm[p^{l-{ 


-u)] 


-{b+/3uj)exp{-ibl) 


cos (Pj^-u) sin(p^M) 


(cj2 - 52 _^ iabUj) p^ 




f sin(p<^[) + 2ipcos{pJ) 





with (3 = duj/v,,, and / = //Aj. Alternatively, the phase was obtained by expansion with respect 

HQ. 



to eigenmodes 
given by 



In particular, the boundary values which determine outside irradiation are 



i 


(6 — Plj) cos{pujl) — i/3 sin(pi^/) 
1 L 


- {b + pto) exp (^-zW j 


J)2 _ ^2 _|_ 

iexp (^—ibfj 


(6 — f3uj) exp 

1 


(1 + P') 
[^bl) - ( 


sin(j9<^/) - 


f 2iPcos{pJ) 

cos ^Pa;^j — "^/^ 


sin (^pj^ 


i~,2 _ ^2 _|_ ^0/^(2; 




(1+^2 


^ sin(p^[) + 2i/3 cos{pJ) 





These cumbersome formulas can be significantly simplified in the regime b^ u, weak dissipation 
at, aq <^ 1, and large impedance mismatch, \P\ <^ 1. In this case, keeping only Fiske resonance 
terms, we obtain 



^^(0) ^ [cos{pJ) - exp{-tbl)]/{ubD), 
e^{-l) ^ [1 - exp(-2W) cos{pJ)]/{CjbD), 
D ^ sin(u;0 + 2(2/3 + al/2) cos(u;/), 



(33) 



In these approximate results we also neglected the term Im [/3(co')], which only slightly shifts res- 
onance positions. Radiation to the right and left, Vl'^^iCj, b), is determined by the values |^a;(0)p 
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and \9^{—l)\'^. At low dissipation, «/ ^ 1, we derive 



^lujulw l-2cos(6/) cos(u;/) + cos2(u;/)±p 



rad - 647^3^2^2 |^|2 

p = 2{al + 2p)[l - cos{bl) cos{ujl)]/b. 



(34) 



The radiation reaches maxima for almost standing Swihart waves at frequencies u = Un = T^n/l 
with n = 1, 2, . . . when cos(6/) cos(tZ'„/) ^ 1. The resonance width is determined by both, the 
dissipation, al, and by the radiation, (3. The perturbation theory is valid in resonance for 16*^ (0) | ~ 
{haluj)^^ < 1 and the radiation power in this linear regime is quite small, Vra.d/w < 10^^ and 
< 1 ;uW/cm at = 10 GHz and 1 THz, respectively. The asymmetry of radiation described by p is 
small. 

Next we derive the dc current at voltage V = $o^/(27rc) and estimate Q. The current / via 
JJ is given by the tunneling quasiparticle contribution. It = atVlw/d, and the Josephson current 
contribution, Ig, 



Is{uj) ^ JJwij 

ij = I W du{cos{ojT — bu)9{T,u))r- 



(35) 
(36) 



The exact result for the reduced Josephson current, ij, is given by 

af,tD/2 



^J 



Im 



(37) 





cos{bl)—cos{p^l) 




{b'^ +pl)sm{pj)-2p^bsm{bl) 


ipl - b^) {uj^ - 62 + iahu) 


[(- 


\- /J^j sin(p^/) + 2iPcos{pJ) 





In the case of weak dissipation the total Josephson current can be split into the dissipation and 
radiation parts. Is = Is, dis + Is, rad- The radiation part plays the same role as dissipation because 
in both cases energy is transfered from the moving vortex lattice to other degrees of freedom (to 
photons in the case of radiation). In the lowest order in A/Aj ^ 1 at the resonance frequencies 
we get 



aSiUJpl 



dis 



rad 



^Qcwal sm.^{bl/2) 



(38) 



2£dcud''""" 327r2AAj62c:;|D|2 
Losses due to radiation are equivalent to those caused by a resistor with Rw = 27r/ {sdU)) attached 
parallel to JJ {Rw ~ 90 ohm-cm for = 1 and u = GHz). The power fed into JJ is P = IV . 
Part of it, [It + Is, dis)^> is dissipated inside JJ, while another part, Ig^ rad^> is radiated. Neglecting 
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non-resonant part, It, we obtain for the radiated fraction, Q = Vrad/V, 

Q = Re[p] — — — (39) 

21 atUJ + I J 

In the regime of weak dissipation and near resonances we approximately obtain 

Q ~ , r = -f = -. (40) 

To clarify the physical meaning of the parameter r, we can represent it as r = Qz{w ^ k~^)J\f 
where 

Qziw » kz') = — ^ « 1 

is the transmission coefficient of the electromagnetic wave at the junction edge (impedance mis- 
match) into free space in the limit wk^ ^ 1 and J\f=l/alis the number of reflections before the 
Swihart wave decays inside JJ. At wk^j ~ 1 and our result is larger than that given by Eq. dU) by 
the factor J\f. As dissipation inside JJ decreases (at and drop), J\f increases. As one can see 
from Eq. (l40b . the relation between the dissipative and radiative dampings is mainly determined 
by competition between the two small parameters, a and d/l. In the linear regime, balu > 1, we 
getA/'< buj. Due to limitation 6 < Aj/A we obtain Q < ((i/A)u)^ ~ 10~^{uj/ujp)'^. 

Figure |3l shows three-dimensional plots illustrating dependencies of the total reduced current 
density J/ Jo = + ij and the radiated fraction Q on the Josephson frequency cu (voltage) 
and field b. For illustration purposes we used toy parameters, I = A, at = 0.01, ag = 0, and 
Po = 0.001. FigurelH shows voltage dependencies of J and Q for several fields. The current shows 
well-known sharp peaks at the Fiske-resonance frequencies, and peak amplitudes oscillate with 
the magnetic field Is^ls^Q]- The conversion coefficient Q shows smooth oscillating behavior 
reaching maxima given by estimate (l40b at the resonances. We also observe another nontrivial 
feature, Q sharply drops at the position corresponding to the Eck resonance conditions, lj = b, 
even though no feature is seen in I-V dependencies for this frequency. For comparison we also 
plotted in Fig. |5] the same dependencies for the case of higher dissipation, at = 0.1. This case 
is quantitatively described by the linear approximation and allows us to trace switching between 
different Fiske peaks with increasing current. 

For a JJ with closed and open edges we obtain 

^^(0) ^-sm{pJ)/{ubD), 

D ^ cos{ul) + (z/2)(/3 + al) sm{ul), (41) 
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FIG. 3: Three-dimensional plots of the current density (left plot) and the radiated fraction (right plot) for 
I = 4, at = 0.01, aq = 0, and /?o = 0.001 obtained within linear approximation. 



1.5 2 2.5 3 




FIG. 4: The voltage dependencies of J and Q at several fields for the same parameters as in the previous 
plot. The current shows sharp increases near the marked Fiske-resonance frequencies. The linear approxi- 
mation does not describe narrow regions near the resonances. Q reaches maxima given by Eq. (l40b at the 
resonances. In the plots for 6 = S.Svr// w 2.75 and Att/I w 3.53 one can see sharp drops of Q at the 
frequencies, corresponding to the Eck resonance condition oj = b. 
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FIG. 5: The voltage dependencies of current J and radiated fraction, Q, at several fields for the higher dis- 
sipation Of = 0.1 (other parameters are the same as in the previous plots). Both J and Q show maxima near 
the mai^ked resonance frequencies. Aitow marks frequency con^esponding to the Eck resonance condition 
uj = b. With cun^ent sweep the system will switch between different peaks, as it is illustrated by the dashed 
lines. 

The Poynting vector of radiation from open edge at x = is given by the expression 

The I-V characteristics defined by Eq. (l35t consists of maxima at the frequencies The part 
of it, where dl/dV < 0, is unstable in the current-biased regime. Experimentally only voltages 
corresponding to the resonance frequencies ujn were observed j^lsl. Numerical calculations 
show that such a behavior occurs in the nonlinear regime, i.e. at very low dissipation. 



liases 



Conclusions for single-junction case 

We have shown that in the case of weakly dissipative junctions Q may become of order unity in 
the linear regime. Nevertheless, the radiation power per unit width, Vra.d/w, given by Eq. (E^ . is 
always small in this regime because the condition < 1 also restricts the power fed into JJ. The 
open question is whether it is possible to get |6'^ | 3> 1 and larger Vra.d in strongly nonlinear regime 
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when dissipation is very low and many Swihart modes are involved. The I-V characteristics in this 
limit have the form of sharp steps (FIske steps) according to experimental data fl^] and numerical 
calculations [Q,!?!. However, the amplitude (pujiO), which determines the radiation power, in highly 
nonlinear regime was not calculated yet. 

In conclusion, accounting for the radiation into the dielectric outside media, we derived the 
dynamic boundary conditions for JJ with the width w much larger than the wave length of the 
electromagnetic wave radiated into the free space. The method of derivation may be extended 
to the case wk^^ <^ 1, and following qualitative conclusions are valid for this case as well. We 
have shown that in the linear regime of Josephson oscillations the radiated fraction Q of the power 
fed into the junction is determined by the number of multiple reflections (i.e. by dissipation rate 
inside JJ) and by the transmission coefficient, Qz, from JJ into free space. Even if Q reaches 
unity, radiation power per unit width of JJ remains small in the linear regime. To probe radiation 
power from JJ in the nonlinear regime, numerical study based on Eqs. (l24b and (l25t is needed. 
We think that measurements of radiation in the best junctions available now, like those studied in 
Ref.H at low temperatures and at intermediate magnetic fields may show higher radiation than 
that observed previously. 



RADIATION IN LAYERED SUPERCONDUCTORS 

For layered superconductors we need to formulate equations for the phase differences (/?„(t, x) 
inside each intrinsic junction n between layers n and + 1 as well as the boundary conditions 
for these variables. Here coordinates of layers are y = {n + l/2)s and layers are parallel to 
the plane {x,z). In general y9„ depend also on the coordinate z, the direction of the applied 
magnetic field Hq, but we will neglect this dependence (straight- vortex approach) assuming that 
width of the crystal along the 2-axis is much larger than the radiation wavelength, k^w ^ 1. The 
derivation of equations for (fnit, x) is similar to that for a single JJ. Using the Maxwell 

equations, the expression for the intralayer supercurrents via the phase of the superconducting 
order parameter inside layers and the Josephson relation for interlayer current one need to exclude 
all variables describing intralayer currents and electromagnetic fields induced by these currents 
by expressing them via Formulation of the boundary conditions is also similar to a single- 
junction case and is based on Eq. (fT^ . However, solution of the equations for phase differences is 
now much more complicated because the vortex structure is two-dimensional (along x and y axis) 
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and may vary depending on the parameters such as the applied dc magnetic field Hq, the length I 
along the x-axis and also on the transport current (velocity of moving lattice). To have significant 
radiation power motion of vortex lattice in different intrinsic junction should be synchronized, as 
in multiple artificial JJ. For that interaction of intrinsic junctions should be strong enough and it 
should favor in-phase vortices in all junctions, i.e rectangular vortex lattice. First, we present the 
equations and the boundary conditions for x) and then discuss solutions for these equations 
and corresponding radiation power in linear regime of Josephson oscillations when perturbation 
theory may be used to solve equations for the phase differences. We will focus on crystals with 
large number of layers on the scale of the London penetration length \ab for intralayer currents. 
In this case super-radiation from many layers becomes possible when there are many junctions 
on the scale of radiation wavelength and when vortices in many junctions are synchronized. We 
will show that in this case significant part of the energy fed into the crystal is converted into the 
radiation. 



Equations for the phase differences 

The crystal length / is along the a;- axis, while the transport current is perpendicular to the layers. 
Thus, the Josephson vortex lattice moves along the x axis, see Fig.0 

As we assumed that there is no dependence of the phase difference, current density j, and the 
electromagnetic fields on z coordinate, the components jz, E-^, B^, and By vanish. The phase 
difference (/5„(t, x) between the layers n and n + 1 is defined as 

27r 

ipn{t,x) = - $„+i - — / dyAy{x), (43) 



ns 



$0 

where is the phase of the superconducting order parameter inside the layer. Intrinsic junction 
n is in the space (n — l/2)s < y < {n + l/2)s. If all intrinsic Josephson junctions are identical 
then the dynamics of the system can be described by reduced coupled equations for the phase 
differences, (pniLx), and reduced magnetic fields hn = By{y = ns)/ Be with Be = $o/(2vrAafeAc) 



(see. e.g., Ref. 114 ill ) 



^ 2 + ""-^ + sinv.„ - — ^ = 0, (44) 
ar^ or au 

' 1 + Va,^ ) + (l + Va,^ ^ = 0. (45) 
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FIG. 6: Schematic picture of layered structure of superconductor and Josephson vortex lattice. The unit 
cell of Bi2Sr2CaCu208 compound and the relation of the crystal structure to the layered model are shown 
in the left. The directions of the applied dc magnetic field Hq, of the dc transport cuiTcnt j, and of the 
radiation Poynting vector Px are also shown. The Josephson vortices forming triangular lattice are shown 
schematically by elliptic cylinders and vi is the velocity of moving lattice. 

Here we used reduced x coordinate, u = x/Xj normalized to the Josephson length Aj = 7s and 
reduced time, r = cupt, where cUp = cj (Acy^) is the plasma frequency, ec is the c-axis high- 
frequency dielectric constant inside the superconductor, \ah and Ac are the London penetration 
lengths, i = \ab/s, s is the interlayer distance, and notates the discrete second derivative op- 
erator, V^v4„ = An+i + v4„„i — 2 An- The dissipation parameters, Vab = ^T^f^ab/ {I'^^c^p) and 
Vc = ^T^f^d i'^c^p), are determined to the quasiparticle conductivities, aab and a^, along and per- 
pendicular to the layers, respectively. The conductivity aab plays the same role as the conductivity 
cTg for a single JJ. Applying the operator — (l + i^ab^) to the first equation and excluding 
hn, we can also derive equations containing only v'„ (m, r) 

with Tc = d/dr+Uc- Eqs. ( B^ represent just charge conservation laws. The first term describes the 
change of electron charges inside the layers and Cooper-pair interlayer tunneling currents, while 
the second term describes superconducting intralayer currents. The terms with the coefficient Uab 
describe quasiparticle dissipative in-plane currents induced by moving Josephson vortices. The 
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static interaction of junctions is described by the term sin while their dynamic interaction is 
described by the term Tc(9/9r)V^v^„. Both are short-range (nearest-neighbor) weak interactions 
and they are not very effective in keeping long-range ordering of vortex lattice along the y-axis. 

In the system of finite number of layers N the dc current with the density J is injected in layer 
1 and extracted from layer N. Then equations for the first and the last junction are obtained from 
Eq. (l46b by putting (po = (pN+i = in linear terms and replacing siny^o = siny^jv+i = j. In a 
finite-layer system the edge junctions differ strongly from other junctions because they have only 
one neighboring junction. 

The electric field inside the superconductor between layers n and n + 1 in terms of the phase 
difference is given as 

^...n+l-^^. (47) 
y,n,n+l ^^^^ V J 

The average (over time or space) electric field determines the Josephson frequency uj. The av- 
erage magnetic field we denote by B and we introduce dimensionless average magnetic field 

b = 2txs\jB/%. 

The parameters of BSCCO at low temperatures are ujp/ (27r) = 0.15 THz, the Josephson critical 
current = 1700 A/cm^, 7 = 500, e, = 12, s = 15.6 A, a^iO) = 2 ■ 10"^ (ohm- cm)-\ 
aab{0) = 4-10'' (ohm- cm)~^ 145] and so £ ^ 100, Uab ~ 0.2 and ^ 5- 10"^. An important feature 
of the high-temperature superconductors is higher relative in-plane dissipation in comparison with 
c axis dissipation, ^ ^^c- Another layered high-T^ compound Tl2Ba2CaCu208 has lower 
anisotropy 7 ^ 150 and as a consequence higher critical current Jc ~ 3 - 10"^ A/cm^ and higher 
plasma frequency, Up/ (27r) ^ 0.75 THz 

Boundary conditions and the Poynting vector of radiation 

We assume that there are only outcoming waves and use Eq. (fT6b for an open edge. In the 
junction n we approximate Bz{y) ~ Bz{yn) and Ey{y) ^ Ey{yn), where yn = sn. This gives the 
following boundary condition at the boundary x = 0, 

^2(0;, n) = f/(u;, n — m)_Ej^(u;, m), (48) 

m 

U{uj,n) ^ (1/2) s[\k^\Jo{k^sn) + ik^No{k^sn)], 7^ 0, 
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FIG. 7: The stack of intrinsic Josephson junctions at x < 0, dielectric at < x < L ^ k^^ and a 
superconducting screen with the London penetration length A., and the thickness D ^ Xg. Such a screen 
introduces additional interaction of vortices in different intrinsic junctions. 

while for n = we need to substitute — (2/7r) ln[l/(|/i;^|s)] for No{kojSn). Finally, we use the 
relations 



Br 



du 



Ey(t,n) 



Br 



J x=0 



Or 



Br 



$0 



J x=0 



2'7lXabXc 



(49) 



to formulate the boundary condition for oscillating part of This gives the following 

general boundary conditions for the phase differences 



s\n - m\) + ik^No{k^s\n - ■m\)](pr, 



du UJec 

In the time representation we obtain again Eq. Wh . as for a single junction. 
For the Poynting vector we derive 

"° 2_,'^o(^<.^s|n-m|)v9„(cj,0)v9;^(a;,0). 



(50) 



64vr3c2Ars 



(51) 



Let us consider now the layered superconductor with a superconducting screen at the right side 
separated from the layered crystal by the dielectric with thickness L small in comparison with the 
wavelength of the electromagnetic wave in this dielectric, see Fig. IT] 

Then, as in the case of a single JJ, we neglect effect of the dielectric on the boundary conditions. 
The solution inside the superconducting screen is given by Eq. (l26t . From this solution we obtain 



dy'K, 



y-y 



Eyiuj,y') 



(52) 
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This relation leads to the boundary condition for the oscillating part of the phase difference 

s(n — m] 



X 



(53) 



Here space derivative of the phase differences and the phase differences are at the point x —0. 
The space derivative now is not small, the parameter {sXj/Xj.) may be of order unity or larger. 
The term in the free energy, J^hc, corresponding to the boundary condition Eq. (l53t . is 



^bc 

w 



dy / dx 



Xs,ujS 



167r3A3 



in 
s{n — m) 



B!{sn] 

4:71 



(54) 



0X(a7, 0), H{an) = Ko{ak)Ko[a{k + n)]. 



This free energy can be compared with the bulk inductive interaction, corresponding to Eq. (l46b . 
seeRef.H 



ind 



duW uVn{u)S/ uVmiu) exp 



\n — m\s 



X 



•ah 



(55) 



Both interactions favor triangular lattice at least in the static case. The coefficient in Eq. (I5?b is 
much larger than that in the inductive free energy, i.e. superconducting screen enhances strongly 
the tendency to form triangular lattice along the y-axis for not very large junction length I. 



Solutions for the phase difference and the radiation power at high fields in large- case 

We consider here the simplest case of layered crystals with large number of layers N > i, 
when we can neglect edge effects along the y-axis, i.e. the difference between the edge (the first 
and the last) and inner junctions. This corresponds to the thickness larger than 0.2 /im with the 
total number of junctions N > 100. 

In the linear regime of Josephson oscillations the general solution for the phase difference has 
the form 

(Pn{r,u) =ujT -bu + K,ri + On{T,u), 6'„(r, u) < 1. (56) 

For Kn = the lattice is rectangular, </?„(r, m) are n-independent. For k„ = nn the lattice is 
triangular, vortices in neighboring layers are in anti-phase positions. For a static lattice such a 
configuration minimizes the energy of the magnetic field inside the crystal with large /, A^. How- 
ever, at small / boundary conditions are inconsistent with triangular lattice at some values of bl. 



and in this case rectangular lattice becomes more favorable ll43h . For moving vortex lattice energy 
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consideration does not work, and here the parameters k„ should be determined by the condition 
that total current / via each junction is the same ll4li . 

JcXjw J du{smLpn{T,u))r = I, ■ (57) 

where Jc = c^o/ (Stt'^sXI) is the Josephson current density. Slowly moving lattice preserves its 
triangular structure. This was confirmed experimentally via observation of magnetic oscillations 
of the flux-flow resistivity with the period of one flux quantum per two junctions Theoretical 



analysis ll4lLl42ll shows that the triangular lattice becomes unstable at the lattice velocity slightly 
smaller than the Swihart velocity cs = cs / {2\ab\^) ■ This instability corresponds to experimen- 
tally observed end point of the first flux-flow branch. It was also shown that interaction with top 
and bottom surfaces leads to significant lattice deformations ll4li . 

Situation at high velocity vi ^ is less clear. Lattice structures and their stability in this 



regime in the case of large lateral size / have been addressed in Ref. 114111 and have been recon 



sidered in Ref. 114211 with the conclusion that for parameters typical for BSCCO no stable regular 
lattice exists at high velocities. Structures and stability of steady states for small lateral sizes / is 
an open issue. 

In the following sections we will estimate the radiation power for rectangular and triangular 
lattice and for lattice with random values of k„. Those are most probable realizations of vortex 
configurations in the large limit. They give also an estimate for maximum radiation power 
which one can anticipate. 

Rectangular vortex lattice 

For rectangular lattice sketched at Fig.[8l = UT—bu+9{T, u), the equation for the oscillating 
part, 9{u), is given by 

For a junction opened at both edges the solution is similar to that of a single JJ, Eq. (l32t . However, 
there is important difference due to the presence of large coefficient £^ in front of the second 
space derivative. This coefficient is due to uniformity of rectangular lattice and corresponding 
supercurrents along the y-axis, see Fig.lU This uniformity leads to large energy of the magnetic 
field, i.e. inductive coupling, Eq. (1551) . resulting in small amplitude of phase variations. To find 
radiation power and I-V characteristics we use results of the perturbation theory for a single JJ. To 
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FIG. 8: Rectangular lattice of Josephson vortices in layered superconductors. The screening currents flow 
along the y axis leading to the large energy of the magnetic field (inductive coupling). Correspondingly, 
the amplitude of phase modulation is weak, but vortices come to the junction edges coherently inducing 
super-radiation. 

make Eq. (l58t similar to that for a single JJ we introduce u = u/£,b = M and / should be replaced 
by l/i. We limit ourself to the frequencies much smaller than the Eck resonance u <^ M. In this 
case the solution has the form 9{t, u) = Re [9^ exp(— icur)] with 

doj{u) ^ e*^"/ (zfe^) + ai exp{ip^u) + 02 exp(-zp^u), (59) 

The coefficients ai and 02 should be found from the boundary conditions given by Eq. ( l5Dt . For 
rectangular lattice at the edge x = we obtain 

= ^sJ2[\kJMk.sn) + tk^NoiKsn)] 9^0). (60) 
Large number of junctions N > {kujs)~^ 

We consider first the case of very large number of layers N > (A;^s)~^. Then 

89^(0) iCj 



du 



^.(0). (61) 



This corresponds io (3 ~ 1/ y^. Important point is that for such (3 the perturbation theory is 
correct for uj\b(!. — Cj\ ^ 1, i.e. practically at all interesting fields and frequencies. We limit 
ourself to the frequencies smaller than the Eck resonance Cj <^h. Thus the amplitude of Swihart 
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waves is proportional to small factor and oscillating part of the phase differences at Nsk^ > 1 
is described by Eqs. (l33t with the extra factor for 9^{0) and 6i^{—l). Further, we need to 
replace ujI by ujl/i ^ 0, while in D we need to put f3 = and a = Vc- Assuming that the 

length of junctions is small enough, ujI <^ i and using results Eq. (l33l) we see that oscillatory 
dependence of the radiation power and of the dc current on ujI drops out and we obtain for the dc 
current at CjM ^ 1 and ujl <^ i the expression without resonances 



I = — — 7 = UcUJ H = — . (62) 

Jcwl hHl OJ 

For the radiation power of N layers, using Eq. (BTT) . we obtain 

Comparing this result with that for a single junction, Eq. (O^t . we see that the additional factor 
N/ (k^si'^) is present. In BSCCO and TBCCO crystals at the frequency = 1 THz the parameter 
kojsi"^ ~ 1. Thus the radiation power of N layers for a moving rectangular lattice is N/{ki^s) 
times larger than from a single intrinsic junction and additional large factor l/{k^s) is due to 
super-radiation in the case Nki^s > 1. 

Let us discuss now the I-V characteristics given by Eq. (l62b . At a given current i > imin we 
have stable solution with positive slope, dl/dV > 0, at a) > cDmin- Here 

/ \l/2 

^niiii = 7^7^ |sin(6//2)|, Vin = Sz/c^^min- (64) 

Hence, moving rectangular lattice cannot exist at currents i < due to super-radiation as power 
fed into the crystal should support radiation as well as quasiparticle dissipation needed to ensure 
stability of this dynamic state. This condition is necessary but it may be not sufficient for the 
stability of moving rectangular lattice. The ratio of radiation power to that of dissipation one is 
r = (tZ'min/'^)^- Thus maximum value of Q is 0.5 at i = imin, and Q drops as current increases 
beyond imin. 

For BSCCO crystals atl = n and 6 = 1 we estimate Umin ~ 4, while Vin ~ 4 x 10"^. At the 
frequency v = 1 THz we obtain V^ad/w ~ N (/iW/cm) for Tl2Ba2CaCu208 crystals. Using the 
energy conservation law, IV = Prad^ we see that to reach this power the total dc current / via the 
junction should be about 0.004 of the critical current. The radiation power in BSCCO crystals at 
the same frequency is about 25 times weaker due to smaller Josephson plasma frequency. 
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Moderate number of junction i < N <^ (k^s) ^ 

Consider now practically more interesting case of moderate number of layers, i < N <^ 
{ki^s)~^. In this case from Eq. (I^UI) we obtain the following boundary condition 

^ = ^«.(0), (65) 



sN 

Pn 



1,1 . 2fc[^ , 
\k,.,\ — I m ■ 



TT I ku, I sN 



This condition is very similar to the boundary condition (l22b for a single junction. 
Solving for oi and 02, we obtain for the edge phase 

cos{ojl/i) — exp{—ibl) 



(66) 



(67) 



sm{ul/i) + i{2(3N + iycl/2i) cos{Col/£) 

Due to Pn <^ 1 at small dissipation we obtain the resonant-type I-V characteristics with Fiske 
resonances atul/i = rm, as in a single JJ. Moreover, for parameters of BSCCO the c-axis quasi- 
particle dissipation is negligible in comparison with the radiation damping. At resonances we 
derive 

1 - (-l)"cos(6n 
ij = z/ecD + — — (68) 

The linear approximation brakes down at 2ibijRePN ~ 1> i-C-j it is valid for number of layers 
N > ^/e^/{b'^CuXab\kuJ\)■ At the frequency ITHz linear approximation is valid if > 100. The 
total radiation power per unit width at the resonance frequency does not depend on number of 
layers and is given by 



rad _ "=c^O"^p 



;-l)"cos(6/) 



(69) 



The independence of this result on A^ appears as a result of cancellation of the A^^ factor in the 
total power of A^ coherently radiating junctions and the factor (3^^ oc A^~^ which determines the 
resonance damping. Away from resonances and for resonances limited by quasiparticle damping 
"Prad oc N^. For BSCCO crystals at a;/27r = 1 THz, we obtain Prad/w^ ~ 24 mW/cm, while for 
Tl2Ba2CaCu208 crystals we estimate "Prad/w ^ 0.5 W/cm. Now tDmin and imin become larger as 
A^ drops: 

-3/2 f ^^/^ 2|sin(6/72)| . 3 

^min = TTT^ , , «min = -^c^xm^- (70) 



''""^ \bHlvJ .jz^;^c ' 2 

As A^ increases beyond N (. the total radiation power does not depend on A^ in Fiske resonances 
and increases oc A^^ outside resonances, while cDmin decrease as N~^^'^. This behavior holds until 
A^ reaches {k^s)^^. After that P^ad increases linearly with A^, while u)min remains constant. 
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It is interesting to compare our estimate with recent large-scale simulations of THz radiation 
out of BSCCO mesa by Tachiki et al. \q\ . They solve the Maxwell equations coupled with the 
equations for the intralayer phases inside the crystal and the dielectric media. They found that 
close to one of the Fiske resonances lattice is mostly disordered but with pronounced rectangular 
correlations promoted by generated standing electromagnetic wave. They found quite powerful 
outside radiation in this state, with power density up to = 3000 W/cm^. Our estimate following 
from Eq. (l69b for 100 junctions occurs to be only slightly smaller ~ 1500 W/cm^. 



Triangular lattice 

Triangular lattice, again in the limit 6^1, behaves quite differently in comparison with rect- 
angular one. The solution has the form 

y9„(r, u) =UJT -bu + 7rn + (-l)"0(r, u) + 6'(r, u). (71) 

Here (p describes the amplitude of the phase oscillations with Josephson frequency u, while 9 



describes oscillations with the frequency 2u; which are induced 



5y moving vortex lattice due to 



nonlinearity of coupled sine-Gordon equations ( B^ . see Ref. |26l According to theoretical esti- 
mates li^ when the lattice velocity approaches the Swihart velocity, the lattice may generate a 
very powerful electromagnetic wave inside superconductor, with power density up to 20W/cm^. 
Unfortunately, this main harmonic of electromagnetic wave at the frequency u experiences full in- 
ternal reflection at the boundary and does not radiate outside. The triangular lattice produce only 
weak outside radiation at frequency 2uj, caused by the homogeneous in c direction component of 
the phase, 9{t,u). In the case of semiinfinite superconductor the radiation power at 2u has been 
estimated in by Artemenko and Remizov [26]. Here we estimate this power for finite-size samples 
when pronounced Fiske resonances are present. 

Putting the solution Eq. (TtTT) into Eq. (l46b we obtain coupled equations for 0(r, u) and 9(t, u): 

+ \ r^TT + T^TT = ~ sm(u;r - bu), (72) 



dr"^ dr \ 4 du"^ J 4 du"^ 
d^9 09 „od^9 ^ ,^ 

— + Uc—-i;— = -(pcos{ujT-bu), (73) 

From the second equation we see that 9 is of order and, as a consequence, terms linear in 9 
were omitted in the first equation in comparison with the term describing dissipation. 
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FIG. 9: Triangular lattice of Josephson vortices in layered superconductors. Screening currents compensate 
each other inside the lattice vortex cell (2s, 1/6) leading to significant reduction of the magnetic energy and 
enhancement of the amplitude of the phase variations in comparison with rectangular lattice. The in-plane 
screening currents lead to intralayer quasiparticles contribution to the dissipation of the moving lattice. Due 
to nonlinearity there are field components at the frequency 2u! with small amplitude. They cause super- 
radiation at the frequency 2uj. 

The boundary conditions for the alternating and homogeneous phase are given by 

m.O, ?m^^,,,eiO). (74) 

for u = 0,1 with P = In the framework of the perturbation theory, at ubatJ ^ 1, the 

amplitude (p, in the lowest order in 1/6, is given by the expression 

4:[cos{q{u + I)) - exp{-ibl) cos{qu)] _ ^~ , 
oq sm[ql) 

where atr = Putting it into the equation for 9 and accounting for boundary conditions 

we obtain 

6 ^ ^[1 - eM-2^bl)] - ' -M-^blW - eM-M)] _ 

2Cu(3£b^^ uj^f3ib^[sm{2ul)+tatJcos{2ul)] 

Keeping only the resonance term, we obtain the radiation power at the frequency 2uj in the linear 
regime 

^c$o^p [cos2(2cD/)-2cos(2cD/)cos(6/'+l][l-cos(M)] 
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Here, as for rectangular lattice, we have large factor {k^s)~^ due to super-radiation at the fre- 
quency 2uj. This differs from the radiation power for rectangular lattice at the same frequency by 
the factor of order unity at u; ~ 1 and 6 ~ 1 in the linear regime. Thus both rectangular and trian- 
gular lattice give approximately the same radiation power for ~ 1 and in both cases the radiation 
is coherent. However, their I-V characteristics and the conversion coefficients Q are very differ- 
ent. The fraction of the power converted into radiation depends on quasiparticle conductivities 
aab and ctc. These conductivities remain nonzero in cuprate superconductors even if temperature 
approaches zero because cuprates are gapless superconductors. As a result, only small part, of 
order (6^£^atrO ^ converted into the radiation in the linear regime. Thus for triangular lattice 
we estimate very small Q oc This translates into strong heating and practically triangular 
lattice cannot provide a realistic source of radiation. 

The I-V characteristics is determined by the behavior of (p and so it is resonant at low dis- 
sipation, as in a single JJ. In particular, for (i) weak dissipation, g ^ 2u) + iat^, and (ii) small 
frequencies u ^ 6/2 we obtain 

cos I 2a)Z ) 



2«b 



-abuj + 



cos{2ujI) — cos{bl) 



sin (2a;/) + (atJ) cos^ ( 2iul 



(78) 



where ai, = Uc + i^abb'^ and atr = ^'c + ^ab^"^- The amplitude of the Fiske peak at 2u;.„/ = nn is 



1 - cos{bl) 



(79) 



Note that triangular lattice excites resonances for the antiphase modes whose frequencies are by 
factor 2£ = 2\ab/s ~ 300 smaller then the frequencies of the homogeneous modes excited by the 
rectangular lattice. Very pronounced Fiske steps have been observed recently by Kim et al. ll23h 

We also estimated that the lattice with random parameters k„ gives incoherent radiation with 
approximately the same power and the conversion coefficient as triangular lattice. 



Conclusions 



We have estimated the radiation power from most probable vortex lattices in the large-iV limit. 
Our estimate for rectangular lattice in BSCCO agrees satisfactory with results of numerical study 
by Tachiki et al. We see that only rectangular lattice gives significant radiation power and 

quite high conversion coefficient, i.e. up to half power fed into the crystal may be converted into 
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the radiation. Thus, not very high currents are needed to get significant radiation power. However, 
the main question to be answered is what are parameters of the stability of rectangular lattice in 
finite-size samples. Other open questions include: 

1. What is the radiation power in highly nonlinear regime for triangular lattice and for rectan- 
gular lattice with moderate number of layers ? 

2. What is the structure of steady states and the radiation power of crystals with not very large 
N < i, when the boundaries along the y-axis become important? 

3. What are boundary conditions and the radiation power in crystals with wk^ ^ 1 ? 

4. Are there fully gapped layered superconducting materials or artificial multilayer system 
which can be used as radiation sources? Presence of finite gap in the whole Fermi surface 
drastically reduces quasiparticle dissipation at low temperatures and the problem of high 
dissipation for the moving triangular lattice in the high-temperature superconductors will be 
avoided. 

5. Is it possible to get stronger radiation from the triangular lattice by modulating the radiated 
edge of the crystal? For example, one can use some periodic layer of dielectric or metal 
between the crystal and the air which may result in more effective conversion of high ky 
electromagnetic fields inside the crystal into low ky outside electromagnetic waves. For 
that the surface modulation has to have the periodicity 2s, i.e., every second layer has to be 
closed. For efficient conversion, one can anticipate that the radiation power would be larger 
by the factor f ^ 10"^ in BSCCO and TBCCO crystals. 
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